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Summary 

During the first reporting period research has concentrated on finishing the modeling 
work required for a representative model of a scramjet propulsion system for hypersonic 
vehicles. An existing hypersonic propulsion code has been adjusted to the Winged-Cone 
Configuration used in this study. In this process the complete force and moment calcula- 
tion has been revised. The advantageous feature of the code to account for angle of attack 
variations was then used to compute the thrust, lift and pitching moment contributions 
of the propulsion system not only for various Mach numbers and fuel equivalence ratios, 
but also for different angles of attack. 
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1 Introduction 


The objective of this research is to address the issues associated with the design of ro- 
bust integrated flight control systems for future hypersonic vehicles with airbreathing 
propulsion systems. It is anticipated that such vehicles will exhibit significant interac- 
tions between rigid body (airframe) dynamics, structural dynamics and engine dynamics. 
The uncertainty in the initial dynamic models developed for these vehicles will also be 
high. The main reason that highly interactive uncertain dynamics are to be expected is 
that scramjet engines will be the primary source of propulsion at hypervelocity speeds. 
Wind-tunnel testing as a result will be limited, and it will be necessary to gain experi- 
ence in actual flight testing of such vehicles. This means that initial flight control system 
design efforts will rely more heavily on theoretical and computer based models, than has 
been the case for subsonic and supersonic aircraft. Also, propulsion system sensitivity to 
angle of attack variations will lead to highly interactive dynamics. 


In this study, the current major research issues from a flight control viewpoint are: (1) the 
development of models that are representative of the interactive dynamics that can occur 
in such vehicles, (2) the development of representative uncertainty models for these dy- 
namics and (3) the development of practical approaches to designing multivariable flight 
control systems that guarantee adequate performance in the presence of uncertainty. The 
research done during the first reporting period has been focussing on item (1). 


The hypersonic vehicle model described in [1] neglects the effects of angle of attack vari- 
ations on propulsion system performance. For long, slender bodies with a considerable 
amount of the compression of the flow going through the engine taking place on the 
forebody, this assumption may not necessarily be valid for the entire flight regime. Fly- 
ing under an angle of attack, an axisymmetric vehicle with engine modules ’wrapped' all 
around the body experiences greater compression on the underside than on the upperside. 
This results in changes in thrust vector magnitude and direction and therefore in flight 
behavior of the vehicle. 


To investigate the angle of attack sensitivity of the propulsion system and the impact 
on vehicle performance the HYTHRUST code described in [2] is used. The code as it 
is presently employed evaluates the thrust, lift, and pitching moment contribution of a 
propulsion unit for various Mach numbers, equivalence ratios, and angles of attack. The 
thrust specific impulse is also computed. 
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Current Research Results 


2.1 Vehicle Description 

As mentioned above, the Winged- Cone Configuration described in [1] is used in this study. 
Main characteristics of this vehicle are an axisymmetric conical forebody, a cylindrical en- 
gine nacelle section with engine modules all around the body, and a cone frustrum engine 
nozzle section. 


The size of an individual engine module was based on information from an input file for 
the ’SRGULL’ program (courtesy of J.D. McMinn, NASA LaRC) and have been modified 
in order to meet certain requirements. These modifications consisted of reducing the 
wedge angle in the inlet in order to avoid subsonic flow inside the scramjet engine at low 
Mach numbers, and reducing the height of an engine module to prevent the bowshock 
from entering the inlet at high Mach numbers and angles of attack. Also, the outer 
shape was somewhat simplified to achieve a better implementation into the code (e.g. 
the streamtubes as they are currently defined do not account for expansion in a vertical 
direction, so this kind of shape was avoided). According to the given width of the inlet a 
total of 32 modules are attached to the vehicle. One of the modules is shown in fig. 1 . 


forebod y aftbody 




Figure 1: Engine module as used in HYTHRUST, all measurements in ft. 
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2.2 Determination of Forces and Moments 

The calculation of the forces and moments created by the propulsion system has been 
revised entirely compared to the version given in [2]. In what follows, a brief outline is 
given describing how the relevant forces and moments are computed. 


The control volume consists of (see fig. 2) 


1. inlet entrance plane 

2. engine surface at body 

3. engine exit plane 

4. engine surface at cowl 

5. aftbody surface. 





Figure 2: Control volume for force computation, vectors indicate outward normal. 

The coordinate system chosen is a two-dimensional windframe system with the x-axis 
positive rearward, the y-axis positive downward, and the pitching moment positive nose 
up. This convention was chosen because the vehicle geometry is already given in this 
system. It is distinguished between two different types of forces: forces created by mo- 
mentum flux, F m , and forces created by pressure, both resolved in directions of the 
coordinates. The pitching moments caused by these forces are M m and A/ p , it is again 
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distinguished if they are caused by an x- or y-force ( M xm , M ym , M xp , M yv ). 


1. Inlet Entrance Plane 

Since for this configuration the streamtubes are parallel to the body x-axis the outward 
normal is fi\ = (— cos a, — sin a), where a is the angle of attack. The incoming flow from 
the forebody has the velocity v\ = i>i(cos 6\, sin 9{). Using the momentum equation and 
the expression 


v-idAi = V\ 77] dA\ 

= — V\ cos(^i — a) dA\ (1) 

the forces caused by the momentum flux for each streamtube are 


dF m \ = pviipidAx) 

= — i?i(cos sin 0i ) dmi 


where 


dm\ — pv\cos{0\ — a)dAi. 
The forces caused by pressure are 

-4 

dF p i = PidAi 

= — pi (cos a, sina)tL4i. 


( 2 ) 

( 3 ) 


( 4 ) 


These force increments have to be summed over all streamtubes. Additionally, they create 
moments around a moment reference center which is assumed for this calculation to be 
at the origin of the coordinate system at the tip of the vehicle. The moments are 


dM = r x dF 

= (x, y) x (dF Xy dF y ) 

= (0,0 ,xdF y — ydF x ) (5) 


So the moments caused by x- and y-forces are 


dM x — -y dF x 
dM y = xdF y , 
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( 6 ) 

( 7 ) 


respectively. Using the expressions above, the resulting moment increments are 


dM xm i 

= yv 1 cos drri\ 

(8) 

dM ym\ 

= — xvisinOidmi 

(9) 

dM xp i 

= ypi cos a dA\ 

(10) 

dA/ ypl 

= — xp\ sin a dA\. 

(11) 


2. Engine Surface at Body 

Since there is no flow across this border of the control volume, only pressure forces are 
relevant. The outward normal is n 2 = (sin a, — cos a), so the forces acting on this surface 
are 


dF p2 = p 2 (sin a, - cos a) dA 2 . (12) 

Therefore, the moments are computed as 

dM xp 2 = — yp 2 sin a dA 2 (13) 

dM VP 2 = — xp 2 cos a dA 2 - (14) 


S. Engine Exit Plane 

For the engine exhaust, the outcoming flow is naturally parallel to the streamtubes and 
therefore to the outward normal n 3 = (cos a, sin a) which leads to the simple expression 

v 3 dA 3 = v 3 dA 3 . (15) 


The forces caused by the momentum flux out of the engine are 

— * 

dF m 2 = V3 ( cos a , sin a) dm 3 


(16) 


where 


dm 3 = pv 3 dA 3 . 


(17) 


The forces caused by pressure are 

— * 

dF p 3 = p 3 (cos a, sin a) dA 3 . (IS) 

The moments can be calculated similarly to the ones at the inlet plane using equations 
(6) and (7) and the expressions for the forces given above. 
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4 . Engine Surface at Cowl 

Like the surface at the body, only pressure forces are relevant. Using the outward normal 
n 4 = (— sin q, cos a), the forces and moments are 


dF p 4 = p 4 (— sin o, cos a) dA 4 (19) 

dM p 4 = p 4 (y sin a, x cos a) dA 4 . (20) 


5. Aftbody Surface 

Again, the aftbody is only affected by pressure forces. Denoting the aftbody angle 
with /? measured in the windframe coordinate system, the outward normal is ns = 
(— sin j 3 , — cos /?) and the forces and moments are 


dF p 5 = ps(- sin /?, - cos /3)dA b (21) 

dMps = p 5 (t/sin/?, -x cos /3)dA 5 . (22) 


Total Forces and Moments 

When summing the forces and moments over the complete engine module it has to be 
considered that the difference in momentum and pressure between the inlet and engine 
exit plane results in a force acting on the fluid. Therefore, the reaction force has to be 
used when forces acting on the vehicle are to be determined. The pressure forces on the 
other surfaces act directly on the vehicle. The total forces are then given by 

Ftot — —(Fmi + F p i + F m 3 + F p3 ) + F P 2 + Fp 4 + F p 5 . (23) 

Likewise, the total moments are 

Mu* = ~{M m i + M p i T M m3 4- Mpf) + M p 2 + M p ± -f M p s- (24) 

So far, the moments were calculated corresponding to a moment reference center at the 
origin of the given coordinate system. For the following calculation of the propulsion 
system database they were transfered to the moment reference center given in [1], 


2.3 Propulsion System Database 

The HYTHRUST code computes propulsion data for one engine module. The Winged- 
Cone Configuration utilizes 32 modules around the whole body. For an axisymmetric 
body at zero angle of attack all modules work under the same conditions and yield an 
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equal amount of thrust, lift, and pitching moment. This changes when the vehicle is flying 
at non-zero angles of attack. Compression will be greatest on the underside and gradu- 
ally decrease when going around the sidewalls to the top. Crossflow components will also 
influence the flow field considerably and alter the flow entering the sidewall engine inlets. 

Since the HYTHRUST code uses an axisymmetric two-dimensional Method of Charac- 
teristics procedure on the forebody, this crossflow cannot be tracked and no answer can 
be given on how the forces will vary for the different engine modules on the sidewalls. 
Implementing a three-dimensional method for the forebody would have been too time- 
consuming so an attempt to estimate the force variation along the sidewalls based on 
experimental results given in [3] and [4] was made. The windtunnel tests described in 
these papers yielded information about the circumferential pressure distribution around 
cones at angles of attack at supersonic speeds. Data for the pressure coefficient cp at the 
the surface for azimuthal angles of 0° (top), 90° and 180° (bottom) was given and an effort 
was made to reproduce these data with the HYTHRUST code. The pressure coefficients 
for top and bottom showed good agreement. However, for an azimuthal angle of 90° our 
modeling approach resulted in exactly the opposite trend from that given in the papers. 
The modeling consisted of running the code with free stream conditions = v ^ cos a 
and introducing an assumed crossflow component v^sina at the 90° inlet. Beyond this 
problem, additional work would have been required to relate the pressure distribution on 
the forebody to the forces and moments created by the propulsion units. 

Since no further information about a reasonably good yet simple estimation could be 
found a linear approach was chosen. The engine data for the whole vehicle are computed 
by evaluating the uppermost and lowermost module, since negligible crossflow will occur 
there, and linearly interpolating for the engine modules in between. 


With the propulsion system forces and moments of the vehicle determined by this proce- 
dure, the thrust coefficient as defined in [1] is 


thrust 


Cp = 


[ft* 


(25) 


where the dynamic pressure is 


2 ' 

1 




The lift coefficient is 


cl = 


lift 

QooSref 




(26) 

(27) 

(28) 
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and the pitching moment coefficient 


cm 


pitching moment 

Qoo^ref C 


H - 


(29) 


The reference area S Te j and reference length c are given in [1]. The fuel specific impulse 
I sp is averaged between the two values for the top and bottom engine module. 


The propulsion system database is computed for flight conditions at an altitude of 100000 
ft (US Standard Atmosphere 1976) and the following conditions: 


Mach numbers: M = 2, 3, 4, 6, 8, 10, 15, 20, 25 

fuel equivalence ratios: <j> = 0.25, 0.5, 0.75, 1.0, 1.5, 2.0, 2.5, 3.0, 4.0, 5.0 

angles of attack: a = 0°, 1°, 2°, 3°, 4° . 

The angle of attack variation is restricted to the condition that there has to be a compres- 
sion on the forebody for the HYTHRUST code to be applicable. Since the forebody cone 
half angle is 5°, the ’ramp angle’ for an angle of attack of o = 4° is only 1° for the flow 
on the uppermost engine module which translates into a very weak compression. Since 
the Winged-Cone Configuration is axisymmetric, the values for negative angles of attack 
correspond to the positive values where lift and pitching moment reverse sign. 


For the zero angle of attack-case data for the thrust coefficient cj and the fuel specific 
impulse I sp already exist [1], Here, an entirely different approach was used to get these 
data, hence they could not be reproduced by the HYTHRUST code. To overcome this 
problem the ’percentage change’ is measured when deviating from a = 0° and this change 
is applied on the existing Winged-Cone data for a = 0° to get the angle of attack varia- 
tion. Using HYTHRUST data, the factor 

ctMo) ~ ctjM = 0 °) f30 x 

ct,h{ q = °°) 

is determined and applied on the Winged-Cone data: 


k T (a) 


ct,wc {<* ) = (1 + fcr(»)) ct,wc{<* = 0°) (31) 

The fuel specific impulse is adjusted similarly whereas lift and pitching moment coeffi- 
cients are not changed. 
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2.4 Results 


To get an impression of how close the HYTHRUST code comes in reproducing the data 
of the Winged-Cone Configuration, a comparison of the thrust coefficients and specific 
impulses for zero angle of attack and four different Mach numbers representing the whole 
flight regime is shown in figs. 3, 4, 5 and 6. Phi is the fuel equivalence ratio, cj is in ft 2 
and I sp is in seconds. 

The HYTHRUST data for M = 2 are considerably lower than the Winged-Cone data 
which can be explained by the fact that in HYTHRUST a pure scramjet engine is mod- 
eled. This type of engine generally does not become effective until the speed increases to 
Mach numbers above five. This trend can be seen for M — 8 in fig. 4 and for M = 15 in 
fig. 5 where the curves for c j and I sp agree more closely with the Langley Winged-Cone 
model. For M = 25 (fig. 6) HYTHRUST yields higher thrust than the Winged-Cone 
configuration which is due to the fact that in the HYTHRUST code no boundary layer 
effects or propulsive drag are considered. 


The subsequent figures 7-26 give an idea of how vehicle performance is affected by the 
propulsion system when the angle of attack is changed. For determining the coefficients, 
the following convention is used: thrust is positive forward and lift is positive upward, 
both expressed in a windframe system. The pitching moment is positive nose up. As 
mentioned above, for cj and I ap the variations due to the angle of attack changes are 
appplied on the given Winged-Cone data for a = 0° in accordance with equations 30 and 
31, whereas the lift and pitching moment coefficients are as given by HYTHRUST. The 
complete propulsion system database for all Mach numbers, fuel equivalence ratios, and 
angles of attack is given in the appendix. 


Thrust and fuel specific impulse are nearly insensitive to changes in angles of attack. Fol- 
lower Mach numbers both cj and I sp slightly decrease for increasing angle of attack, and 
for higher Mach numbers they may decrease as well as increase marginally depending on 
the fuel equivalence ratio used. 


The lift added to the vehicle by the propulsion system increases for increasing angle of 
attack. The magnitude of the lift coefficient first increases, then decreases again when the 
Mach number is varied from 2 to 25. It can be seen especially for M = 15 (fig. 20) that 
there is a certain ’kink’ or even a ’jump’ for M = 25 (fig. 24) in the c^-curve for small 
equivalence ratios. This is the result of a simplified modeling of the combustion process 
as a ’simple’ Rayleigh line heat addition. It is assumed that when reaching a certain 
maximum temperature in the combustion chamber no more fuel is burned and thrust is 
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Mach 

cm q f° r basic vehicle 

c\f a for prop, system 



<£ = 0.5 

<£=1 

4 

-0.00175 

-0.000119 

-0.000091 

6 

-0.00100 

-0.000451 

-0.000334 

10 

0.00050 

-0.000550 

-0.000515 

15 

0.00100 

-0.000535 

-0.000401 

20 

0.00110 

-0.000342 

-0.000063 

25 

0.00110 

-0.000174 

0.000213 


Table 1: c\j a comparison for or = 1° 

increased by injecting excessive hydrogen. The transition into this fuel-rich environment 
causes the described effect which is also visible in the graphs of the pitching moment. 


The greatest impact on vehicle performance by the propulsion system is in the pitch- 
ing moment. For small Mach numbers the pitching moment contribution is throughout 
positive (nose up) (see fig. 9) and increases linearly with fuel equivalence ratio. With 
increasing speed, cm becomes negative first for small, and then for all equivalence ratios 
(see figs. 13 and 17). This trend is reversed when the Mach number is increased further 
(fig. 21) and for M — 25 cm is positive again for almost all equivalence ratios (fig. 25). 
Furthermore, the contribution of the propulsion system to the vehicle pitching moment 
due to angle of attack variations shows a stabilizing tendency. Table 1 gives a comparison 
of the pitching moment sensitivity to angle of attack changes 


c M a 


d£M_ 

da 


(32) 


for the basic vehicle (derived from fig. 21 in [l] ) and for the propulsion system for the 
equivalence ratios <£ = 0.5 and <£ = 1 evaluated at a = 1°. It can be seen that for Mach 
numbers around 10 the influence of the propulsion system reaches the same magnitude as 
the basic vehicle’s influence but is reversed in sign and therefore has a stabilizing influence. 
This effect starts at around M=6 and can still be observed at M=15. 


2.5 Conclusions 

The investigation shows that the impact of propulsive variations due to angle of attack 
changes on vehicle performance is not negligible for typical hypersonic vehicles. In this 
case, it is not so much that thrust is affected but that the pitching moment displays a 
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significant sensitivity for this type of long and slender vehicle. Even small variations in 
the forebody conditions can result in considerable changes in the pitching moment which 
make the vehicle behavior highly unpredictable. This indicates the presence of uncer- 
tainty for further control studies. 


3 Future Research 

During the next reporting period research will concentrate on developing the linear mod- 
els of the aircraft using the created propulsion sytem model. In this effort we will rely on 
the results provided by NASA using the POST program as well as investigate the possi- 
bility of employing a different code performing the same or a similar task. Furthermore, 
the uncertainty modeling will be carried out and a paper for the 1993 AIAA Guidance, 
Navigation and Controls conference will be prepared that will illustrate a control system 
design based on a robust performance specification. 
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Figure 3: Comparison HYTHRUST and Winged Cone data, M = 2 
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Figure 8: Lift variation, M = 2 
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Thrust, M ■ 4 



Figure 11: Thrust variation, M = 4 

xlO- 3 Lift, M «=4 



Figure 12: Lift variation, M = 4 




















Figure 16: Lift variation, M = 8 















Figure 20: Lift variation, M = 15 
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Figure 24: Lift variation, M = 25 
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Appendix 

A.l Input Data File HYTH.DAT 


The input data file used by the HYTHRUST code specifying freestream properties as well 
as geometrical and combustion data is given below. The expressions on the right hand 
side of the file represent the names of the variables the values of which are given on the 
left hand side. 


8. 

408.6 23.198 0. 

AMinf ' 

Tinf Pinf alpha 

1.4 

1716. 

Gamma ] 

Rgas 


3 


NPF 



0 . 

0 . 0 . 

XF(1) 

YF(1) 

CF(1) 

73.0 

6.39 .15658 

XF (2) 

YF (2) 

CF (2) 

150. 

13.122 .0762 

XF (NPF) 

YF (NPF) 

CF (NPF) 

FALSE 


Lout 



2.536 

4.5 4.5 

H XW1 

XW2 


3. -3 

. -3. 3. 

D1 D1ST D2 

D2ST 

150. 

14.95 159. 14.95 

XCL YCL 

XEX 

YEX 

5. 

0 0 . 

XCC NSTRUT ALST 

1.38 

1716. 1.2E9 .96 .0292 

G2 RG2 

QR 

EtaB FCRIT 

.5 0 

. 

ER AM 



3 


NPA 



159. 

13.122 0.0762 

XA(1) 

YA(1) 

CA(1) 

180. 

9.796 0.10206 

X A(2) 

YA(2) 

CA(2) 

200. 

6.63 0.15082 

XA(NPA) 

YA(NPA) 

CA(NPA) 
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A. 2 Propulsion database 


Here, the complete database for the propulsion system’s thrust, lift, and pitching moment 
contribution including fuel specific impulse data is given for all combinations of Mach 
number, fuel eqivalence ratio and angle of attack. The expressions in the first block rep- 
resent the given variables and explain the structure of the database. 


M 

equ. 

ratio 



alpha 

ct 

cl 

cm 

Isp 

alpha 

alpha 

ct 

cl 

cm 

Isp 





2.0 

0.25 




0.0 

205.042 

0 . 000E+00 

0 . 000E+00 

5672.40 

1.0 

204.627 

0 . 116E-02 

0 . 640E-04 

5668.14 

2.0 

203.439 

0 . 232E-02 

0 . 124E-03 

5658.67 

3.0 

201.568 

0 . 347E-02 

0 . 176E-03 

5644.48 

4.0 

198.800 

0 . 467E-02 

0 . 186E-03 

5623.08 

2.0 

0.50 




0.0 

315.572 

0 . 000E+00 

0 . 000E+00 

4365.00 

1.0 

314.950 

0. 116E-02 

0 . 767E-04 

4361.78 

2.0 

313.169 

0.232E-02 

0 . 149E-03 

4354.61 

3.0 

310.363 

0 . 348E-02 

0 . 214E-03 

4343 . 85 

4.0 

306.212 

0.469E-02 

0 . 237E-03 

4327.63 

2.0 

0.75 




0.0 

391.292 

0 . 000E+00 

0 . 000E+00 

3608.30 

1.0 

390.541 

0.117E-02 

0 . 893E-04 

3605.68 

2.0 

388.387 

0 . 233E-02 

0 . 175E-03 

3599.85 

3.0 

384.996 

0 . 350E-02 

0 . 252E-03 

3591.07 

4.0 

379.977 

0.471E-02 

0.288E-03 

3577 . 87 

2.0 

1.00 




0.0 

451.446 

0 . 000E+00 

0 . 000E+00 

3122.20 

1.0 

450.602 

0.117E-02 

0 . 102E-03 

3119.97 

2.0 

448.176 

0 . 234E-02 

0.200E-03 

3115.00 

3.0 

444.361 

0 . 351E-02 

0 . 290E-03 

3107.51 

4.0 

438.710 

0.472E-02 

0 . 338E-03 

3096.25 
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2.0 

1.50 




0.0 

553.891 

0 . 000E+00 

0. OOOE+OO 

2553.90 

1.0 

552.905 

0 . 118E-02 

0 . 127E-03 

2552.14 

2.0 

550.066 

0 . 236E-02 

0 . 251E-03 

2548.18 

3.0 

545.607 

0 . 353E-02 

0.366E-03 

2542.21 

4.0 

538.995 

0.475E-02 

0 . 439E-03 

2533.25 

2.0 

2.00 




0.0 

651.470 

0 . OOOE+OO 

0. OOOE+OO 

2252.90 

1.0 

650.363 

0 . 119E-02 

0 . 153E-03 

2251.39 

2.0 

647.172 

0 . 237E-02 

0 . 301E-03 

2248.00 

3.0 

642.164 

0.356E-02 

0.442E-03 

2242.85 

4.0 

634.733 

0.479E-02 

0.541E-03 

2235.15 

2.0 

2.50 




0.0 

753.045 

0. 000E+00 

0. OOOE+OO 

2083.20 

1.0 

751.823 

0. 120E-02 

0 . 178E-03 

2081.84 

2.0 

748.289 

0 . 239E-02 

0 . 352E-03 

2078.79 

3.0 

742.752 

0 . 358E-02 

0 . 518E-03 

2074.13 

4.0 

734.526 

0.482E-02 

0 . 642E-03 

2067.18 

2.0 

3.00 




0.0 

861.741 

0. OOOE+OO 

0. OOOE+OO 

1986.60 

1.0 

860.402 

0.120E-02 

0 . 203E-03 

1985.34 

2.0 

856.522 

0 . 240E-02 

0.403E-03 

1982.49 

3.0 

850.450 

0 . 360E-02 

0 . 594E-03 

1978.15 

4.0 

841.421 

0.485E-02 

0.743E-03 

1971.67 

2.0 

4.00 




0.0 

1104.097 

0. OOOE+OO 

0. OOOE+OO 

1909.00 

1.0 

1102.517 

0 . 122E-02 

0 . 254E-03 

1907.85 

2.0 

1097.920 

0 . 244E-02 

0 . 504E-03 

1905.23 

3.0 

1090.741 

0 . 365E-02 

0.746E-03 

1901.21 

4.0 

1080.045 

0.491E-02 

0.946E-03 

1895.25 

2.0 

5.00 




0.0 
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1860.77 

4.0 

111.892 

0 . 103E-01 

-0.221E-02 

1865.98 

10.0 

1.50 




0.0 

130.855 

0. OOOE+OO 

-0 . 935E-09 

1404.00 

1.0 

130.635 

0 . 253E-02 

-0 . 454E-03 

1404.41 

2.0 

129.913 

0 . 507E-02 

-0 . 917E-03 

1405.37 

3.0 

128.848 

0.782E-02 

-0 . 148E-02 

1406.37 

4.0 

128.020 

0 . 104E-01 

-0 . 197E-02 

1409.82 

10.0 

2.00 




0.0 

140.495 

0. OOOE+OO 

0. OOOE+OO 

1130.60 

1.0 

140.323 

0 . 255E-02 

-0.399E-03 

1130.89 

2.0 

139.742 

0 . 509E-02 

-0 . 805E-03 

1131.54 

3.0 

138.931 

0 . 786E-02 

-0 . 131E-02 

1132.17 

4.0 

138.490 

0 . 104E-01 

-0 . 174E-02 

1134.61 

10.0 

2.50 




0.0 

147.949 

0. OOOE+OO 

0. OOOE+OO 

952.45 

1.0 

147.827 

0 . 256E-02 

-0.344E-03 

952.66 

2.0 

147.394 

0 . 512E-02 

-0 . 694E-03 

953.12 
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3.0 

146 . 844 

0 . 790E-02 

-0 . 114E-02 

953.53 

4.0 

146.790 

0 . 105E-01 

-0 . 150E-02 

955.34 

10.0 

3.00 




0.0 

154.732 

0 . 000E+00 

0. OOOE+OO 

830.11 

1.0 

154.658 

0 . 257E-02 

-0 . 289E-03 

830 . 27 

2.0 

154.369 

0 . 515E-02 

-0 . 582E-03 

830.60 

3.0 

154.075 

0 . 794E-02 

-0 . 967E-03 

830.87 

4.0 

154.395 

0 . 106E-01 

-0 . 127E-02 

832.25 

10.0 

4.00 




0.0 

168.765 

0 . 000E+00 

0. OOOE+OO 

679 . 04 

1.0 

168.784 

0 . 260E-02 

-0.179E-03 

679.13 

2.0 

168.768 

0 . 520E-02 

-0 . 360E-03 

679.31 

3.0 

168.958 

0 . 802E-02 

-0.625E-03 

679.40 

4.0 

169.996 

0 . 107E-01 

-0.796E-03 

680.27 

10.0 

5.00 




0.0 

184.954 

0 . OOOE+OO 

0. OOOE+OO 

595.35 

1.0 

185.062 

0 . 263E-02 

-0 . 696E-04 

595.41 

2.0 

185.309 

0 . 526E-02 

-0 . 137E-03 

595.50 

3.0 

185.968 

0 . 811E-02 

-0 . 284E-03 

595.49 

4.0 

187.714 

0 . 108E-01 

-0.327E-03 

596.06 

15.0 

0.25 




0.0 

28.204 

0. 000E+00 

0. OOOE+OO 

1948.70 

1.0 

28.191 

0. 133E-02 

-0.392E-03 

1947.96 

2.0 

28.129 

0 . 267E-02 

-0 . 783E-03 

1946.36 

3.0 

27.957 

0 . 401E-02 

-0 . 118E-02 

1944.10 

4.0 

27.562 

0 . 532E-02 

-0 . 157E-02 

1941.70 

15.0 

0.50 




0.0 

50.137 

0. OOOE+OO 

0 . 237E-09 

1732.10 

1.0 

50.078 

0.200E-02 

-0 . 536E-03 

1733.77 

2.0 

49.862 

0 . 398E-02 

-0 . 107E-02 

1738.90 

3.0 

49.381 

0 . 596E-02 

-0. 161E-02 

1746.35 

4.0 

48.495 

0 . 785E-02 

-0.213E-02 

1756.53 

15.0 

0.75 




0.0 

64.495 

0. OOOE+OO 

-0.475E-09 

1485.40 

1.0 

64.369 

0 . 204E-02 

-0 . 447E-03 

1488.26 
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2.0 

63.949 

0.405E-02 

-0 . 888E-03 

1496.75 

3.0 

63.105 

0.601E-02 

-0.132E-02 

1508.92 

4.0 

61.709 

0.783E-02 

-0 . 172E-02 

1524.46 

15.0 

1.00 




0.0 

74.624 

0 . 000E+00 

0 . 000E+00 

1289.00 

1.0 

74.521 

0 . 205E-02 

-0.405E-03 

1291.31 

2.0 

74.163 

0.407E-02 

-0 . 802E-03 

1298.15 

3.0 

73.392 

0 . 604E-02 

-0.119E-02 

1307.94 

4.0 

72 . 043 

0.788E-02 

-0 . 154E-02 

1320.43 

15.0 

1.50 




0.0 

88.168 

0 . 000E+00 

0 . 000E+00 

1015.30 

1.0 

88.132 

0 . 207E-02 

-0 . 320E-03 

1016.89 

2.0 

87.957 

0.410E-02 

-0 . 630E-03 

1021.59 

3.0 

87.452 

0 . 610E-02 

-0 . 928E-03 

1028.29 

4.0 

86.383 

0 . 796E-02 

-0 . 119E-02 

1036.83 

15.0 

2.00 




0.0 

97.358 

0 . 000E+00 

0 . 000E+00 

840 . 85 

1.0 

97.393 

0 . 209E-02 

-0 . 235E-03 

842.01 

2.0 

97.421 

0.414E-02 

-0.458E-03 

845.45 

3.0 

97.220 

0 . 616E-02 

-0.666E-03 

850.34 

4.0 

96.502 

0.804E-02 

-0 . 829E-03 

856.54 

15.0 

2.50 




0.0 

104.769 

0 . 000E+00 

0 . 000E+00 

723.89 

1.0 

104.872 

0 . 210E-02 

-0 . 151E-03 

724.78 

2.0 

105.096 

0.418E-02 

-0 . 287E-03 

727.42 

3.0 

105.195 

0 . 621E-02 

-0 . 403E-03 

731.16 

4.0 

104.830 

0 . 812E-02 

-0.472E-03 

735.90 

15.0 

3.00 




0.0 

111.594 

0 . 000E+00 

0.475E-09 

642.53 

1.0 

111.762 

0 . 212E-02 

-0 . 662E-04 

643 . 24 

2.0 

112.173 

0.422E-02 

-0 . 115E-03 

645.35 

3.0 

112.557 

0 . 627E-02 

-0 . 140E-03 

648.32 

4.0 

112.530 

0 . 820E-02 

-0 . 116E-03 

652.08 

15.0 

4.00 




0.0 

125.637 

0 . 000E+00 

0.356E-09 

542 . 54 
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1.0 

125.930 

0 . 216E-02 

0 . 103E-03 

543.04 

2.0 

126.697 

0.430E-02 

0 . 228E-03 

544.50 

3.0 

127.622 

0 . 639E-02 

0 . 385E-03 

546 . 55 

4.0 

128.230 

0.836E-02 

0 . 596E-03 

549.13 

15.0 

5.00 




0.0 

141.733 

0 . 000E+00 

0 . OOOE+OO 

489 . 64 

1.0 

142.151 

0.220E-02 

0 . 272E-03 

490.02 

2.0 

143.273 

0.437E-02 

0 . 572E-03 

491.13 

3.0 

144.733 

0 . 651E-02 

0.910E-03 

492.67 

4.0 

145.960 

0 . 852E-02 

0 . 131E-02 

494.59 

20.0 

0.25 




0.0 

24 . 548 

0 . 000E+00 

0. OOOE+OO 

1885.10 

1.0 

24.499 

0 . 104E-02 

-0 . 277E-03 

1883.26 

2.0 

24.335 

0 . 197E-02 

-0.514E-03 

1879.99 

3.0 

23.888 

0 . 300E-02 

-0.789E-03 

1872.01 

4.0 

23.046 

0.402E-02 

-0 . 107E-02 

1856.84 

20.0 

0.50 




0.0 

44.459 

0 . 000E+00 

0. OOOE+OO 

1707.00 

1.0 

44.319 

0 . 153E-02 

-0 . 368E-03 

1708.57 

2.0 

43.861 

0 . 293E-02 

-0 . 685E-03 

1714.39 

3.0 

42.810 

0.440E-02 

-0 . 105E-02 

1720.75 

4.0 

41.036 

0 . 574E-02 

-0 . 137E-02 

1725.44 

20.0 

0.75 




0.0 

58.078 

0 . 000E+00 

0. OOOE+OO 

1486.60 

1.0 

57.820 

0 . 129E-02 

-0.146E-03 

1490.83 

2.0 

56.975 

0.243E-02 

-0.239E-03 

1503.27 

3.0 

55.210 

0.367E-02 

-0.396E-03 

1519.21 

4.0 

52.419 

0.482E-02 

-0 . 565E-03 

1535.19 

20.0 

1.00 




0.0 

68.082 

0 . 000E+00 

0 . 134E-09 

1307.00 

1.0 

67.825 

0 . 130E-02 

-0.894E-04 

1310.38 

2.0 

66.967 

0 . 245E-02 

-0 . 126E-03 

1320.32 

3.0 

65.109 

0 . 370E-02 

-0 . 229E-03 

1333.02 

4.0 

62.111 

0.487E-02 

-0.346E-03 

1345.69 


20.0 


1.50 
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0.0 

82.204 

0 . 000E+00 

0 . 334E-09 

1052.10 

1.0 

81.981 

0. 133E-02 

0 . 232E-04 

1054.40 

2.0 

81.195 

0 . 250E-02 

0 . 990E-04 

1061.14 

3.0 

79.353 

0.377E-02 

0 . 106E-03 

1069.70 

4.0 

76.252 

0.497E-02 

0 . 929E-04 

1078.15 

20.0 

2.00 




0.0 

92.411 

0 . 000E+00 

-0 . 134E-09 

887.02 

1.0 

92.234 

0 . 135E-02 

0 . 136E-03 

888.68 

2.0 

91.561 

0 . 255E-02 

0 . 324E-03 

893.57 

3.0 

89.828 

0 . 385E-02 

0.441E-03 

899.74 

4.0 

86.779 

0 . 506E-02 

0 . 532E-03 

905.75 

20.0 

2.50 




0.0 

100.970 

0 . 000E+00 

0 . OOOE+OO 

775.34 

1.0 

100.839 

0 . 138E-02 

0 . 249E-03 

776.61 

2.0 

100.282 

0 . 260E-02 

0 . 549E-03 

780.34 

3.0 

98.677 

0 . 392E-02 

0.776E-03 

785.02 

4.0 

95.717 

0 . 516E-02 

0 . 970E-03 

789.52 

20.0 

3.00 




0.0 

108.965 

0 . 000E+00 

0. 000E+00 

697.28 

1.0 

108.878 

0 . 140E-02 

0 . 361E-03 

698.29 

2.0 

108.432 

0.265E-02 

0 . 774E-03 

701.25 

3.0 

106.947 

0 . 399E-02 

0 . 111E-02 

704.95 

4.0 

104.073 

0 . 526E-02 

0 . 141E-02 

708.46 

20.0 

4.00 




0.0 

125.283 

0 . 000E+00 

-0 . 167E-09 

601.27 

1.0 

125.276 

0.145E-02 

0 . 587E-03 

601.97 

2.0 

125.029 

0 . 274E-02 

0 . 122E-02 

604.02 

3.0 

123.750 

0.414E-02 

0 . 178E-02 

606.54 

4.0 

121.000 

0 . 545E-02 

0 . 229E-02 

608.87 

20.0 

5.00 




0.0 

143.513 

0 . 000E+00 

0. OOOE+OO 

551.01 

1.0 

143.581 

0. 150E-02 

0.812E-03 

551.54 

2.0 

143.516 

0 . 284E-02 

0 . 167E-02 

553.09 

3.0 

142.402 

0.429E-02 

0 . 245E-02 

554.96 

4.0 

139.703 

0 . 565E-02 

0.316E-02 

556.64 
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25.0 

0.25 




0.0 

12.904 

0 . 000E+00 

0. OOOE+OO 

1096.50 

1.0 

12.822 

0 . 738E-03 

-0.140E-03 

1095.23 

2.0 

12.556 

0 . 148E-02 

-0.286E-03 

1091.19 

3.0 

12.083 

0.226E-02 

-0.461E-03 

1082.93 

4.0 

11.415 

0 . 308E-02 

-0 . 672E-03 

1067.45 

25.0 

0.50 




0.0 

24.635 

0 . 000E+00 

-0 . 107E-09 

1046.70 

1.0 

24.451 

0 . 107E-02 

-0 . 170E-03 

1047.05 

2.0 

23.865 

0 . 212E-02 

-0 . 349E-03 

1047.84 

3.0 

22.788 

0 . 285E-02 

-0.413E-03 

1045.53 

4.0 

21.241 

0 . 333E-02 

-0.401E-03 

1036.80 

25.0 

0.75 




0.0 

33.842 

0 . 000E+00 

0. OOOE+OO 

958.60 

1.0 

33.529 

0 . 563E-03 

0 . 189E-03 

961.18 

2.0 

32.540 

0 . 115E-02 

0 . 345E-03 

968.06 

3.0 

30.807 

0 . 181E-02 

0.411E-03 

974.18 

4.0 

28.457 

0.259E-02 

0 . 358E-03 

973.49 

25.0 

1.00 




0.0 

41.566 

0 . 000E+00 

-0.855E-10 

883.00 

1.0 

41.207 

0 . 577E-03 

0 . 256E-03 

885.11 

2.0 

40.066 

0 . 118E-02 

0.476E-03 

890.73 

3.0 

38.116 

0 . 185E-02 

0 . 610E-03 

897.72 

4.0 

35.508 

0 . 257E-02 

0 . 649E-03 

903.26 

25.0 

1.50 




0.0 

54.489 

0 . OOOE+OO 

0 . 192E-09 

771.68 

1.0 

54.067 

0 . 606E-03 

0 . 390E-03 

773.18 

2.0 

52.717 

0 . 123E-02 

0.740E-03 

777.15 

3.0 

50.396 

0 . 193E-02 

0 . 993E-03 

782.04 

4.0 

47.294 

0 . 268E-02 

0 . 114E-02 

785.80 

25.0 

2.00 




0.0 

65.602 

0. 000E+00 

0 . 855E-10 

696.81 

1.0 

65.135 

0 . 635E-03 

0 . 524E-03 

697.94 

2.0 

63.634 

0 . 129E-02 

0 . 100E-02 

700.92 

3.0 

61.042 

0 . 202E-02 

0 . 138E-02 

704.56 

4.0 

57.579 

0 . 279E-02 

0.163E-02 

707.26 
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25.0 

0.0 

1.0 
- 2.0 

3.0 

4.0 

25.0 

0.0 

1.0 

2.0 

3.0 

4.0 

25.0 

0.0 

1.0 

2.0 

3.0 

4.0 

25.0 

0.0 

1.0 

2.0 

3.0 

4.0 


2.50 

75.820 

0 . OOOE+OO 

0. OOOE+OO 

644.28 

75.317 

0 . 663E-03 

0 . 659E-03 

645 . 17 

73.689 

0.135E-02 

0 . 127E-02 

647.51 

70.868 

0 . 210E-02 

0.176E-02 

650.33 

67.099 

0 . 290E-02 

0 . 212E-02 

652.36 

3.00 

85.595 

0. 000E+00 

0. OOOE+OO 

606.11 

85.059 

0.692E-03 

0.793E-03 

606.83 

83.315 

0 . 140E-02 

0 . 153E-02 

608.73 

80.283 

0 . 218E-02 

0 . 214E-02 

610.99 

76.234 

0 . 300E-02 

0.260E-02 

612.55 

4.00 

104.682 

0. OOOE+OO 

0. OOOE+OO 

555.95 

104.083 

0 . 750E-03 

0 . 106E-02 

556.46 

102.116 

0.152E-02 

0 . 206E-02 

557.80 

98.680 

0 . 235E-02 

0 . 290E-02 

559.35 

94.091 

0 . 322E-02 

0 . 358E-02 

560.31 

5.00 

123.848 

0. OOOE+OO 

0. OOOE+OO 

526.19 

123.185 

0.807E-03 

0 . 133E-02 

526.58 

120.994 

0.163E-02 

0 . 258E-02 

527.58 

117.151 

0 . 252E-02 

0 . 367E-02 

528.70 

112.019 

0 . 344E-02 

0.456E-02 

529.31 
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